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Abstract 

Progress in computational atmospheric dynamics is 
exhibiting the ability of numerical simulation to de- 
scribe instability processes associated with turbulence 
observed at altitudes between 15 and 25 km in the lower 
stratosphere. As these numerical simulation tools ma- 
ture, they can be used to extend estimates of atmo- 
spheric perturbations from the present gust database 
for airplane design at altitudes below 15 km to alti- 
tudes between 25 and 50 km where aerospace plane 
operation would be at hypersonic speeds. The amount 
of available gust data and number of temperature per- 
turbation observations are limited at altitudes between 
15 and 25 km. On the other hand, in-situ gust data at 
higher altitudes are virtually nonexistent. The uncer- 
tain potential for future airbreathing hypersonic flight 
research vehicles to encounter strong turbulence at 
higher altitudes could penalize the design of these vehi- 
cles by undue cost or limitations on performance. Be- 
cause the atmospheric structure changes markedly with 
altitude, direct extrapolation of gust magnitudes and 
encounter probabilities to the higher flight altitudes is 
not advisable. This paper presents a brief review of 
turbulence characteristics observed in the lower strato- 
sphere and highlights the progress of computational at- 
mospheric dynamics that may be used to estimate the 
severity of atmospheric transients at higher altitudes. 

Nomenclature 

CAT clear air turbulence 

Fr Froude number, Np r 

g acceleration caused by gravity, m/sec 2 
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HAT high-altitude turbulence 

HICAT high-altitude clear air turbulence 

h geopotential altitude, km 

h rate of change of altitude, m/min 

L height of the flow barrier in 

the Froude number 

MSL mean sea level 

PSD power spectral density 

p pressure, mb 

rms root mean square 

T temperature, K 

U horizontal windspeed component 

Ude derived equivalent gust velocity 

U second derivative of U with respect to 

altitude 

USAF United States Air Force 

w vertical velocity component 

x horizontal coordinate 

8 atmospheric static stability parameter, 

g 66/06h 

6 derivative operator 

A finite difference operator 

0 potential temperature, T(1000/p)° 286 

p atmospheric density 

Introduction 

Advanced hypersonic research vehicles, such as the 
National Aero-Space Plane, require discrete atmo- 
spheric perturbation models for design at altitudes 
above the present gust and turbulence design criteria 
based on in-situ measurements. Recent developments 
indicate that numerical simulation of small-scale at- 
mospheric dynamics could assist in the specification 
of these discrete atmospheric perturbation models for 



hypersonic aircraft design, especially at high flight alti- 
tudes. The atmospheric turbulence environment speci- 
fied for design of civil and military aircraft is described 
as a function of altitude in terms of discrete and ran- 
dom gust inputs. 1 '’ 2 These discrete inputs are given 
as the derived equivalent gust velocity. 3 

Continuous random inputs are specified by a power 
spectral density (PSD) shape, wavelength scale factor, 
and root-mean-square (rms) value. 4 These design data 
are based on extensive gust loads surveys aboard op- 
erational aircraft and on true gust velocity measure- 
ment programs on specially instrumented military and 
research aircraft. 5-8 The adequacy of these criteria 
for gust loads design of subsonic aircraft to altitudes 
of approximately 15 krn has been demonstrated on 
narrow- and wide-body transport aircraft. 5 For the 
higher altitudes, gust loads and true gust velocity data 
have been obtained by U-2 aircraft (Lockheed Corpo- 
ration, Burbank, California) to altitudes of approxi- 
mately 20 km. 9-10 These data have also been ob- 
tained by supersonic research aircraft to altitudes of 
25 km. 11-14 

Because hypersonic flight research will use altitudes 
to approximately 50 km, these existing in-situ atmo- 
spheric turbulence data cover approximately one-half 
of the altitude range required. Moreover, subsonic 
transport airplane design has emphasized gust input 
effects on the limit load and fatigue life design. In 
contrast, hypersonic vehicles are more likely to require 
emphasis on propulsion and flight control systems re- 
sponses. Response of hot structures to high-altitude 
turbulence and gusts will also be required. In flight 
regimes where more than one system is sensitive to 
the perturbation inputs, the interacting responses will 
probably involve more than one of the input pertur- 
bation components. For example, critical responses 
may simultaneously involve vertical and horizontal gust 
components as well as combinations of gust velocity 
and ambient density or temperature. In this regard, a 
chief limitation of the present gust database and design 
criteria is that they have only been established for sin- 
gle, independent perturbation component inputs. That 
is, for present conventional aircraft design, combina- 
tions of vertical and horizontal gusts are prescribed as 
being statistically independent. In addition, present 
aircraft design criteria do not link either the horizon- 
tal or vertical gust components with perturbations in 
temperature, density, or pressure. 

These concerns for hypersonic aircraft design apply 
not only to the present database altitudes below 25 km 
but also become even more critical at altitudes from 
25 to 50 km where higher speed airbreathing aero- 
dynamic flight will be pioneered. Thus, two major 
aspects of the turbulence environment remain to be 
explored for hypersonic aircraft. The first aspect is 


statistical characterization of the intensity and amount 
of turbulence throughout the stratosphere. The second 
aspect involves description of the nature of coupling be- 
tween gust motion components and pressure, density, 
and temperature state perturbations. Formulation of 
higher altitude design criteria and perturbation mod- 
els addressing both of these aspects requires use of ad- 
vanced numerical simulation techniques for small-scale 
atmospheric perturbations. In addition, appropriate 
observations to validate the simulations for these ap- 
plications are needed. 

This paper provides a brief overview of the char- 
acteristics of turbulence observed in the lower strato- 
sphere. Current studies for the use of small-scale, two- 
dimensional, numerical atmospheric dynamics simula- 
tions for specification of the higher altitude perturba- 
tion environment are also discussed. Finally, a simple, 
generic, discrete atmospheric perturbation model pro- 
posed for assessment of hypersonic aircraft sensitivity 
to combinations of atmospheric gusts and thermody- 
namic transients is described. 

Metric units have been used throughout this paper. 
For the convenience of readers who are more familiar 
with the U.S. Customary System, conversion factors 
are as follows: 


U.S. 

Metric 

Conversion factor 

Celsius, °C 

Kelvin, K 

"C = K - 273.15 

Fahrenheit, °F 

Kelvin, K 

K = (5/9) (°F + 459.67) 

Feet, ft 

Meters, m 

m = ft x 0.3048 

Millibar, mb 

Hectopas- 
cal, hPa 

mb = hPa 

Nautical 
miles, n. mi. 

Kilome- 
ters, km 

n. mi. = km x 0.54 
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Atmospheric Structure 

Figure 1 shows the atmospheric layers with respect 
to the average January conditions at Edwards, Cali- 
fornia. These data are based on rawinsonde observa- 
tions to 30-km altitude and nearby rocketsonde data 
for altitudes from 30 to 70 km and serve as a repre- 
sentative example of midlatitude winter conditions. 15 
In the troposphere, temperature generally decreases 
about 6.5 K/km altitude, and windspeed generally 
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increases up to the tropopause level at the base of 
the stratosphere. In the lower stratosphere, tempera- 
ture tends to be nearly constant with altitude (isother- 
mal), and windspeed generally decreases to minimum 
values at altitudes between 20 or 25 km. In the 
upper stratosphere, the rate of temperature increase 
reaches approximately 2.8 K/km. Windspeeds increase 
with altitude through the upper stratosphere into the 
mesosphere, where the temperature again decreases 
at rates of 2 to 4 K/km altitude. The division be- 
tween the stratosphere and mesosphere is termed the 
stratopause. 

In the lower troposphere, sensible heat is actively 
exchanged by mixing processes which bring air into 
contact with the earth and sea surfaces. Low alti- 
tude warming predominates during the day at lower 
latitudes. Conversely, cooling prevails at night and at 
higher latitudes. Much heat is also exchanged through- 
out the troposphere by moisture phase changes. Solar 
radiative heat input to the atmosphere is located pre- 
dominantly in the upper stratosphere where chemical 
species, such as ozone, absorb solar ultraviolet radi- 
ation. Between these two heat source regions, tem- 
peratures are cooler in the upper troposphere and 
lower stratosphere. Basically, all atmospheric layers 
are subject to infrared radiative heat loss which, at the 
surface, may be notably limited by cloud cover. 

Because air is compressible, temperature changes ex- 
perienced by atmospheric airparcels are not limited to 
heat gain or loss. These changes are also affected by 
adiabatic expansion and compression when airparcels 
rise and descend. As airparcels are displaced adiabat- 
ically, that is, without heat gain or loss, up or down 
in the atmosphere, the temperatures change by a rate 
of nearly 10 K/km of altitude displaced. Cooling oc- 
curs for upward displacements, and warming results 
from downward displacements. Thus when vertical mo- 
tion is not uniformly distributed in the atmosphere, 
airparcels with the most relative displacement expe- 
rience the greatest temperature differences from the 
surrounding airmass. Subsequently as rising airparcels 
become cooler and more dense than the surrounding 
airmass, the force of gravity causes buoyancy effects 
to restore the airparcels toward their original altitude 
levels. Conversely, as descending airparcels become 
warmer or less dense than the surrounding ambient air, 
these parcels also experience restoring buoyancy which 
forces them toward their original altitude. 

Because ambient temperatures increase with altitude 
in the stratosphere, the airparcel temperature contrast 
grows at rates of 10 to 13 K/km displaced, causing the 
buoyancy forces to act as a stiff spring. In comparison, 
weaker spring force effects are experienced in the tro- 
posphere and mesosphere, where airparcel temperature 
contrast grows at rates from 4 to 8 K/km displaced. 


This effect, in which airparcel altitude displacements 
are opposed by buoyancy forces, is much stronger in 
the stratosphere than in either the troposphere or the 
mesosphere. As a result, the term stratosphere is quite 
appropriate for this deeply stratified region of the at- 
mosphere. The differences between wind and tempera- 
ture structure in the troposphere and stratosphere also 
lead to inherent differences in atmospheric turbulence 
characteristics. Differences in turbulence characteris- 
tics are believed to include gust magnitudes, length 
scales, horizontal and vertical patch dimensions, patch- 
iness or spacing between patches, and characteristics of 
associated temperature perturbations. As indicated in 
the next section, these characteristics have only been 
partially observed to date. 

Turbulence Observations in the Lower 
Stratosphere 

Early operation of the U-2 airplane in the 1950’s 
provided gust loads data to altitudes near 23 km 
over several representative areas of the world. 9 The 
USAF High-Altitude Clear Air Turbulence (HICAT) 
Program subsequently conducted true gust velocity 
measurements in the stratosphere with the U-2 airplane 
over various terrain and geographical areas. 10 The 
XB-70 (North American Rockwell, Los Angeles, Cal- 
ifornia) and YF-12 (Lockheed Corporation, Burbank, 
California) prototype supersonic cruise aircraft also 
provided high-altitude gust acceleration data and true 
gust velocity measurements as a part of aeronautical 
flight research programs. 11-14 

Results from these flight programs (Fig. 2) gener- 
ally indicate that the amount of turbulence expressed 
as a percentage of total flight distance decreases from 
a maximum near the jetstream and tropopause alti- 
tudes to a minimum near altitudes of 20 or 25 km. 
Figure 2 shows that the fraction of miles in clear air 
turbulence (CAT) for use with PSD criteria also de- 
creases significantly in the lower stratosphere. 7 The 
altitudes at which the amount of turbulence decreases 
most changed somewhat from one sampling program 
to another. Presumably, these changes in the altitude 
of turbulence decrease were caused by seasonal differ- 
ences in the atmospheric structure as well as geograph- 
ical area. 13 

The number of gust loads or normal acceleration 
peaks cumulated from large to small and experienced as 
a function of overall flight distances also decreases dra- 
matically (Fig. 3) at the higher altitudes of the super- 
sonic airplane data. The cross-hatched area on Fig. 3 
shows that the range of gust loads acceleration experi- 
ence for subsonic transport jet aircraft is close to the 
supersonic aircraft experience at altitudes between ap- 
proximately 12 to 17 km. In contrast at altitudes above 
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18 km, gust accelerations experienced in supersonic 
flight are markedly less than those experienced in ei- 
ther subsonic or supersonic flight at lower altitudes. 
Because 25 km is generally near the altitude of minimal 
windspeed, gust loads are expected to remain relatively 
mild at altitudes between 20 to 25 km. At higher alti- 
tudes, however, winds generally increase although the 
direction may reverse. Such changes in windspeed indi- 
cate that the amount and intensity of turbulence may 
be expected to increase at altitudes of 30 to 50 km in 
the upper stratosphere. 

Some insight into the observed decrease in turbu- 
lence at altitudes in the lower stratosphere is provided 
by inspection of seasonal trends, geographical patterns, 
and empirical associations with attending meteorolog- 
ical features. Figure 4 shows the seasonal variation 
in the relative amount of turbulence encountered in 
the lower stratosphere by the YF-12 airplane on flights 
made from Edwards, California. 13 These data are 
smoothed by a 3-month moving mean for the altitude 
layer from 12 to 17 km in the lower stratosphere. They 
depict a maximum in the winter and early spring when 
lower altitude jetstream and airmass frontal activity is 
greatest and a minimum in the summer and early fall 
when the winds in the lower stratosphere are weak and 
have reversed to easterlies. 15 

Meteorological Features of Turbulence 
in the Lower Stratosphere 

On individual days, meteorologists have noted the 
association of high-altitude turbulence (HAT) at alti- 
tudes above 12.5 to 25 km with mountain-wave activ- 
ity and strong winds blowing over the tops of large 
cumulonimbus clouds. 10,16-17 The association with 
wave activity was examined further for several of the 
XB-70 turbulence encounter days. 18-19 These findings 
not only supported the indications from studies by 
HICAT meteorologists 10,17 but also demonstrated the 
association of HAT with both rawinsonde balloon rise 
rates and with a parameter important to mountain- 
wave behavior, that is, the curvature of the lower al- 
titude wind profile. Nominal rawinsonde balloon rise 
rates are typically from 250 to 300 m/min. These nom- 
inal rise rates may change gradually with balloon ma- 
terial characteristics, amount of inflation, and atmo- 
spheric temperature. Variations around the nominal 
rise rate will sometimes occur as a result of up- or 
down-drafts associated with cloud or wave activity in 
the atmosphere. 

Large and cyclic variations in the balloon rise rates 
are indicative of mountain waves and turbulence en- 
countered in the stratosphere. 18 Examples of the bal- 
loon rise rate variations associated with high-altitude 
turbulence conditions are shown by the profiles in 


Figs. 5(a) to (d). 18 These profiles represent two bal- 
loon ascents for a day with mountain- wave conditions 
and HAT and two ascents for another day with nei- 
ther mountain waves nor HAT. For HAT cases, nearby 
balloon rise rate variations in the troposphere are gen- 
erally in the range of 50 to 150 m/min. Nonturbu- 
lent cases generally have rise rate variations from 30 
to 70 m/min. This study encompassed 112 balloon as- 
cents from 7 locations over 39 high-altitude flight days. 
The criteria for large and cyclic balloon rise rate vari- 
ations correctly indicated 80 percent of the HAT en- 
counters on these days. 

Another parameter associated with HAT is an in- 
crease in the wind profile curvature between altitudes 
from 3 to 7.6 km. Curvature is the second derivative 
of windspeed, U /U, with respect to altitude divided by 
windspeed. A curvature increase between two altitudes 
(At/ jU > 0) is conducive to the formation of large am- 
plitude mountain waves. 19 These findings were devel- 
oped by identifying areas of expected mountain- wave 
activity and those with positive curvature parameter 
changes in the lower troposphere between layers cen- 
tered near 3 and 7.6 km and examining these areas 
with respect to turbulence encountered along the flight 
track of the XB-70 airplane (Fig. 6). For evaluation of 
the forecast, skills provided by the curvature parameter 
nonturbulence regions were also specified in smooth 
areas of the flight track where altitude changes were 
minimal. The flight example shown in Fig. 6 has two 
HAT encounters in areas of expected mountain-wave 
activity where the curvature criteria were fulfilled, one 
smooth segment of the track in another area of ex- 
pected wave activity, and four smooth track segments 
which were outside of these areas. Overall evaluation 
of this forecast parameter on 15 high-altitude flights 
resulted in correct identification of 40 turbulence en- 
counters out of a total of 47 encounters and one false 
alarm in 43 nonturbulence areas. Thus, the role of 
lower altitude wave activity has been empirically es- 
tablished by independent investigators for a significant 
portion of HAT cases encountered by both subsonic 
and supersonic aircraft. This finding was significant 
but not surprising because mountain waves had been 
known to increase the severity of CAT at passenger 
jet cruise altitudes. Many researchers have been in- 
strumental in the analysis and demonstration of CAT 
enhancement by mountain- wave- induced vertical dis- 
placement of shear layers to cause Kelvin -Helmholtz 
wave amplification and instability. 20-25 

These convincing demonstrations of the Kelvin- 
Helmholtz instability role in CAT were undoubtedly 
a large part of the answer to explaining the genera- 
tion of many CAT encounters. However, the Kelvin- 
Ilelmholtz instability explanations did not suffice as 
a large part of the solution to the CAT forecast and 
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avoidance problem because its application would rely 
on forecasting the vertical displacement of windshear 
layers. Accurate forecasts of windshear layer displace- 
ments with the necessary detail in time and space are 
not feasible. The related atmospheric processes often 
involve such local phenomena as mountain waves with 
small-scale atmospheric variations which entail subgrid 
phenomena. None the less, mountain- wave field studies 
and theoretical analyses of the attending atmospheric 
dynamics have made considerable progress in the con- 
cepts of mountain- wave propagation behavior. 26-28 

An example of upward propagating mountain-wave 
activity which causes turbulence in the lower strato- 
sphere is depicted by the temperature and wind fields 
analyzed for a case of in-si tu aircraft observations in 
Fig. 7. 28 Moderate and severe turbulence was encoun- 
tered at several aircraft sampling altitudes as denoted 
by the “hats” in the figure. Potential temperature is 
the temperature that the airparcel would have if adi- 
abatically compressed to 1000 mb, that is, approx- 
imately 100 m mean sea level (MSL) pressure alti- 
tude. In Fig. 7, potential temperature contours are 
denoted by solid lines. Horizontal wind contours (iso- 
tachs, lines of constant speed) are shown as dashed 
lines. Upstream windspeeds exceed 60 m/sec near 7.5- 
and 12-km altitudes. At higher altitudes, the upstream 
winds generally decrease to less than 20 m /sec at 20-km 
altitude. The potential temperature contours, which 
approximate trajectory streamlines, show the downs- 
lope winds over the mountain ridge. Successive wave 
crests indicate partially trapped wave components at 5 
to 7 km in the middle troposphere. 

Although wave amplitudes are diminished near 
11 km, a large amplitude jump is generated at the 
14- to 17- km level. At 16-km altitude, upstream of 
the jump, the zero isotach indicates a narrow region of 
flow reversal. Such reversals suggest that cooler, more 
dense air from lower altitudes has been brought above 
warmer air to form a local layer of convective overturn- 
ing. Unfortunately, the measurement scale was not suf- 
ficiently detailed to indicate the cell sizes of instabilities 
involved in this region. Strong horizontal temperature 
gradients are noted in the jump and downstream of 
the peak. At 16 km, the horizontal windspeed increases 
from less than zero in the blocked or overturning region 
to 40 m/sec in the peak of the jump. The windspeed 
then decreases to between 20 and 30 rn/sec in a pro- 
nounced horizontal temperature gradient with strong 
turbulence. Near 19 km, the highest measurement al- 
titude, light turbulence was still observed, and wave 
amplitudes and wavelengths were decreased. Wave 
cases are seldom identical and often change measur- 
ably over periods of an hour or more. Just the same, 
the features of this case are representative of wave 
character changes from the lower troposphere to the 


tropopause and lower stratosphere when moderate and 
severe stratospheric turbulence occurs. Note that the 
sharp density changes resulting from strong temper- 
ature perturbations in similar cases are often on the 
order of 5 percent. Such changes would have measur- 
able impacts on aerodynamic coefficients and propul- 
sion system mass flow for an aerospace plane. 

Numerical Simulation of Wave 
Propagation 

Relevant progress in applying numerical simula- 
tion to atmospheric gravity wave propagation, par- 
ticularly into the lower stratosphere, has been re- 
cently reviewed. 29 Several numerical simulation codes 
are now available and may have a role in advancing 
the quantitative understanding of small-scale perturba- 
tions observed in the lower stratosphere and in the pre- 
diction of disturbances at altitudes above the present 
airplane design criteria database. 30-33 Numerical tech- 
niques are particularly applicable to the study of per- 
turbations in the lower stratosphere because quanti- 
tative treatment of mountain-wave behavior by closed 
form solution is limited to simple analytical profiles 
of wind and thermal structures and to small obsta- 
cles or flow barriers. Multiple layers must be used to 
represent the irregular atmospheric wind and temper- 
ature profiles frequently associated with most cases of 
strong turbulence. Because the juncture between the 
troposphere and stratosphere at the tropopause and 
jetstream level is also an inherent, first-order, atmo- 
spheric discontinuity, multiple layers are also involved 
for estimating the amplitude of waves propagating into 
the stratosphere. Therefore, it is appropriate to seek 
numerical solution techniques to adequately account 
for the influence of irregular temperature and wind pro- 
files on mountain-wave propagation and turbulence in 
the stratosphere. 

An initial experiment to ascertain the viability of 
numerical simulations for application to small-scale 
stratospheric perturbations is described in the next 
three subsections. 34 A brief overview of other related 
numerical simulation studies follows in a later section. 

Intercomparison Cases 

In this initial comparison experiment, observed cases 
of disturbances in the lower stratosphere were com- 
pared with numerical simulation codes having vari- 
ous underlying physical assumptions, numerical formu- 
lations, and topographical representation schemes. 34 
Simulation input consisted of the topographic or bar- 
rier description as the lower boundary and the up- 
stream profiles of atmospheric temperature and wind 
as the input condition at every vertical gridline. This 
comparison experiment was initiated to identify dif- 
ferences between the numerical codes that would be 
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important to applications in the stratosphere and to 
provide a practical assessment of the capability for 
specifying the atmospheric profiles and topography for 
the numerical simulation initial conditions in applied 
cases. Moreover, agreement between observation and 
models indicates the viability of using the numerical 
simulation models to specify perturbation characteris- 
tics at altitudes above available observations. 

Six cases of documented in-situ wave observations in 
the stratosphere 27 ” 28, 35 ~ 37 were selected for the pilot 
project intercomparisons of simulations by the numer- 
ical codes. 30 ” 33 These simulation codes can run with- 
out restrictions of linearity or hydrostatic equilibrium. 
Differences in these codes include treatment of com- 
pressibility, topographic representation, grid structure 
(fixed cartesian versus variable-grid resolution and or- 
thogonal curvilinear), and viscid effects. 

Three topographic reliefs relative to the barrier are 
represented by the following cases: 

• Sierra- Nevada and White Mountains 

• Rocky Mountain Continental Divide west of Den- 
ver, Colorado 

• Rocky Mountain ridgeline east of Alamosa, 
Colorado 

The Sierra-Nevada and White Mountains have lower 
terrain upstream than downstream. The Rocky Moun- 
tain Continental Divide has higher terrain upstream 
than downstream. The Rocky Mountain ridgeline has 
approximately the same terrain elevation upstream as 
it does downstream from the barrier. 

All six cases had in-situ aircraft observations at mul- 
tiple altitudes from the upper troposphere to 18 km in 
the stratosphere. In four cases, strong turbulence was 
encountered in the stratosphere. The other two cases 
exhibited well-established wave activity in the lower 
troposphere but did not report notable turbulence or 
wave activity above 14 km. In the atmospheric layer at 
and below the mountain ridge level altitude, the cases 
represent an inverse square root of the Froude number 
(Np r = U**2/Lg). These cases range from approx- 
imately 0.8 to 3.5 and, therefore, cover more than a 
narrow range of dynamic flow conditions near the to- 
pography. This application to the atmosphere differs 
slightly from the Froude number as conventionally used 
in hydrodynamics. Here, U is the windspeed upstream 
of a mountain ridge height, L, and the atmospheric sta- 
bility ( gj3 or g60/66h) is used in place of g alone because 
of atmospheric buoyancy. 29 The analyzed atmospheric 
wind and temperature data for the wave propagation 
event depicted in Fig. 7 is an example of the in-situ 
data used for comparison with the numerical simula- 
tions for these cases. 28 


Case Preparation for Numerical Model Input 

Upstream atmospheric wind and temperature pro- 
files were prepared for input to the numerical models 
from routinely archived upper air network rawinsonde 
data and synoptic charts. Figures 8 and 9 show ex- 
amples of the original upper air observations and the 
resulting temperature and wind input profiles. For the 
first stage of the intercomparisons, preparation of the 
input atmospheric profiles and topographic representa- 
tion was accomplished independently from the numer- 
ical modeling. These atmospheric data were interpo- 
lated in time and space to give representative values of 
wind and temperature at 1-km altitude intervals. As 
suggested by the scatter of the observed data around 
the input profiles, this data value selection process was 
subjective with a tendency to produce constant wind 
shear and lapse rate layers over large altitude segments. 

Topographic barrier descriptions for the simu- 
lation lower boundary condition were also speci- 
fied by a three-way, subjective compromise among 
aircraft ground track, ridge layer wind direction, 
and topographic irregularities. These independently 
prepared input data were then given to the numerical 
modeling personnel at the Universities of Washington 
(Seattle) and California (Los Angeles) without identifi- 
cation of the dates or in-situ case observations. Imple- 
mentation of these input data in the models and selec- 
tion of grid resolution, time step, model domain, and 
output simulation times were not rigidly restricted. In- 
stead, such decisions were left to the expert judgment 
of the individual modelers. 

Preliminary Comparison Results 

Results of the initial intercomparisons between sim- 
ulations and observed wave and turbulence were highly 
encouraging. 34 As an example, Figs. 10(a) to (d) show 
the strong contrast between simulations produced for 
wave perturbation cases, with and without strong tur- 
bulence in the stratosphere, for the potential tempera- 
ture and horizontal windspeed fields. All three models 
which simulated these two cases exhibited similar con- 
trast between the one with and the one without turbu- 
lence in the stratosphere. This turbulence case simu- 
lation corresponds to the in-situ measurement analysis 
shown in Fig. 7. The simulation replicates the dimin- 
ished wave amplitudes near 11 km, potential temper- 
ature jump between 14 and 17 km, large horizontal 
windspeed changes at 16 km, and decreased wave am- 
plitudes and wavelength above 19 km. 

Preliminary assessment of the simulation results for 
the six cases was based on inspection of qualitative 
wave propagation features. The qualitative features 
included wave character in the troposphere, wave am- 
plitude propagation behavior in the tropopause zone, 
and wave-breaking tendency indicative of turbulence 
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in the lower stratosphere. Inspection of the analyzed 
aircraft observation data for the 6 cases yielded 15 of 
these qualitative features for evaluation of the initial 
outputs available from 2 models. For the 30 qualita- 
tive wave propagation features examined on the simu- 
lations, agreement was judged as good for 21, as fair 
for 4, and as poor for 5 of the evaluated features. 

Quantitative comparison of amplitudes for the larger 
wave perturbations at the flight data altitudes were 
also assessed. Quantitative examination emphasized 
the predominant perturbation wave heights, vertical 
velocity, and windspeed changes. Magnitudes of these 
features were generally estimated from contour plots 
for the observed and the numerical model output data. 
More than 50 percent of the 92 comparisons agreed to 
within a factor of 2. 34 In judging the goodness of this 
agreement, the following factors were considered: 

• Perturbation magnitudes vary by more than a fac- 
tor of 2 between subjective turbulence intensity- 
rating categories. 

• Aircraft measured perturbation magnitudes typi- 
cally change by 10 to 100 percent for repeated data 
runs at given altitudes within a 3-hr period. 

• Manual data extraction from contour presenta- 
tions is not precise. 

• Strong perturbations are not expected to be 
steady state in either the atmosphere or the nu- 
merical simulations. 

In view of these considerations, the present intercom- 
parison results are judged to be excellent and very en- 
couraging for the application of numerical simulation 
to wave perturbation phenomena in the stratosphere. 34 

Further assessment will consider the sensitivity of 
the individual models and comparisons among them to 
the following parameters: 

• Resolution and smoothing of topographic features 
as well as representation by grid blocks (restricted 
to horizontal and vertical surfaces) versus terrain- 
following constructions 

• Computational grid resolution in horizontal and 
vertical space dimensions as well as in the time- 
step size 

• Duration of the simulations in atmospheric time 

• Resource requirements in terms of computer time, 
cost, and user expertise 


These assessments will consider present and future ap- 
plications which may include extension to altitudes 
above 50 km and finer grid resolutions to study 
the influence of shear layers with Kelvin-Helmholtz 
instability. 

Related Wave Propagation Studies 

Generation of atmospheric turbulence in the strato- 
sphere is basically the result of buoyant imbalance and 
motion stresses much the same as in the lower atmo- 
sphere. In the troposphere, however, buoyant imbal- 
ance is frequently caused by heating of the ground or 
clouds, releasing of latent heat in condensation pro- 
cesses, mixing as winds blow over rough terrain, or 
strengthening of upper air wind and temperature gra- 
dients. In contrast, a larger portion of turbulence in the 
stratosphere is caused by wave processes triggered by 
mountains, cloud barriers, or irregular jetstream dy- 
namics which amplify as they propagate upward. In 
the troposphere, vertical propagation of some wave en- 
ergy is frequently prevented by airmass layers. These 
layers have negative or neutral static stability, that is, 
no buoyant restoring force. Propagation is also pre- 
vented by reflection caused by strong changes in wind- 
speed with altitude on the underside of the jetstream 
or by significant directional shifts in the wind profile. 
As wave motion components which escape the afore- 
mentioned limitations propagate upward, amplitude 
and velocity perturbations of these components tend 
to grow approximately with the inverse square root of 
the density ratio between altitudes. 29,38 Thus, upward 
wave leakage amplifies until it generates its own insta- 
bilities or meets other limiting factors in the wind or 
temperature static stability profile. 

Ground-based measurements and satellite observa- 
tions are among methods used to gather evidence 
of wave amplification and instability with propa- 
gation to higher altitudes in the mesosphere and 
thermosphere. 39 ” 41 The extremely low density at these 
altitudes leads to large amplitude disturbances which 
are observable by remote means. These phenomena 
are believed to explain much of the atmospheric den- 
sity variations inferred from space shuttle entry tra- 
jectory data. 42 The most dramatic of these was the 
alleged “density hole” which strongly perturbed the 
Space Shuttle Columbia causing a jolting response sim- 
ilar to heavy turbulence during atmospheric entry near 
76-km altitude on STS-4. Most concern for airbreath- 
ing hypersonic flight research is, however, at much 
lower altitudes where dynamic pressures will be suf- 
ficient for acceleration and level cruise. Thus, the need 
for realistic specification of atmospheric perturbations 
for hypersonic flight is most critical at altitudes from 
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25 to 50 km. These altitudes arc above the limit of 
most in-si tu data, yet they are still somewhat below 
most of the remote observation studies. 

The mechanisms involved in wave propagation from 
the middle stratosphere to the mesosphere have been 
studied by a limited number of numerical simulations 
which have required considerable computer time. 43 
These simulation studies have illustrated the ampli- 
fication and instability processes and have provided 
useful additions to conventional theoretical interpre- 
tations. As an example. Figs. 11(a) to (d) show the 
potential temperature contours resulting from a nu- 
merical simulation of wave propagation between alti- 
tudes of 30 and 50 km. 43 The wave response is trig- 
gered by a barrier located 50 km downstream from the 
inflow boundary and at the base altitude of the sim- 
ulation, 30 km. The potential temperature contours, 
which approximate flow streamlines for nearly steady 
conditions, show the wave amplitude increases with al- 
titude. At 10-hr simulation time (Fig. 11(a)), the up- 
ward motion becomes very steep in patches between 
40- and 50-km altitude. In local regions where the 
potential temperature contours slant back to the left 
as they go up, the simulation produces colder, more 
dense air over warmer, lighter air. The prevailing left- 
to-right wind decelerates and begins to reverse. Con- 
tinuation of this process (Figs. 11(b) and (c)) leads to 
further reversal of the wind in these local regions and 
to wave breaking or overturning as the small areas of 
cold air accelerate downward. Turbulence and mixing 
in these small areas develop regions of nearly constant 
potential temperature as noted by the tendency to form 
pockets of separated potential temperature isotherms 
(Fig. 11(c)) at 14 hr into the simulation. Finally at 
16 hr (Fig. 11(d)), the wave perturbations and turbu- 
lence decay before the onset of another cycle of wave 
propagation, growth, and overturning. 

In spite of the value of this upper stratosphere nu- 
merical simulation work, many questions of wave be- 
havior remain to be answered in the lower stratosphere 
(below 30 km) before the knowledge of higher altitude 
wave behavior can contribute to design risk definitions 
for hypersonic vehicles. These questions involve the 
frequency of significant wave component generation at, 
or propagation through, the tropopause. The critical 
layer behavior in the lower stratosphere must also be 
ascertained. Studies of the instability zone associated 
with the critical level will greatly assist in this effort. 

Wave propagation studies have recently addressed 
the lower stratospheric structure where wind decreases 
with altitude. Theoretically, critical levels occur where 
the windspeed equals the wave-phase speed (or zero 
windspeed for a standing wave). Results of this work 
emphasize the following facets of wave propagation 
behavior: 


• Critical level behavior, wherein propagating wave 
energy is converted to turbulent and kinetic en- 
ergy, is commonly manifested a few kilometers 
below the actual critical level altitude where 
the primary upward propagating wave meets re- 
flected wave components which are propagating 
downward. 44 

• Critical level criteria do not need to be completely 
fulfilled to induce a layer of instability in the lower 
stratosphere. 45 That is, the windspeed decrease 
with altitude does not actually need to reach zero 
(or reverse) with respect to the wave. 

• Because of the effects of Earth rotation, longer 
wavelength inertial-gravity wave propagation can 
induce instabilities at multiple levels over a broad 
or deep altitude zone depending on the atmo- 
spheric profile structure and mode wavelength or 
wave number. 29,45 

• Nonhydrostatic effects impact the momentum flux 
and cause wave perturbations in the stratosphere 
to occur further downwind than predicted by 
solutions based on the simplifying hydrostatic 
assumption. 46 

Figures 12(a) to (d) show an example from the previ- 
ously cited studies. 44 This example depicts instability 
development in a deep layer of decreasing wind with an 
upward propagating wave component approaching the 
critical level. Wave propagation behavior in this sim- 
ulation study corresponds to what may be expected 
in the altitude band from 15 to 30 km at times. In 
this case, the input upstream windspeed decreases with 
constant shear from 10 m/sec at 15 km (the base 
altitude of the layer) to 0 m/sec at 25-km altitude 
(the critical level) and continues in reverse direction 
to —5 m/sec at 30-km altitude. Figures 12(a) to (d) 
show the streamfunction resulting from monotonic ex- 
citation having a wavelength of 10 km and an ampli- 
tude less than 150 m introduced at the base altitude of 
15 km of the nonlinear, spectral simulation. For clarity 
of illustration, Figs. 12(a) to (d) show the streamfunc- 
tion patterns for only the 20- to 26-km altitude section 
of the simulation at times of 4, 9, 10, and 12 hr after 
introduction of the 10-km wave at the lower bound- 
ary. In addition to delineate the instability pattern, the 
streamfunction contour interval resolution is decreased 
by a factor of 10 just below the 22-km altitude. 

The first phase, depicted in Fig. 12(a) at 4 hr into 
the simulation, shows the wave crests and troughs slope 
upstream with altitude as they propagate to higher al- 
titudes. As the propagation process continues, critical 
level interactions reflect wave energy which modifies 
the form of the upward propagating wave. This inter- 
action subsequently produces sloping zones of reverse 
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flow at altitudes between 23 km and the critical level 
at 25 km as shown in Fig. 12(b) for 9 hr into the simu- 
lation. At this time, rapid development of smaller scale 
instabilities is initiated and becomes evident between 
altitudes of 22 and 24 km just 1 hr later as shown in 
Fig. 12(c). As the instability processes enter their later 
stages at 12 hr into the simulation time (Fig. 12(d)), 
note that the small-scale disturbances not only extend 
below 22 km but some also begin to propagate above 
the critical level at 25 km. This instability zone be- 
low the critical level, somewhat above the tropopause 
and jetstream altitudes, will often be at altitudes where 
airbreathing vehicles accelerate from supersonic to hy- 
personic speeds. 

This progress in numerical wave propagation and 
simulation studies indicates that atmospheric dynamics 
can have a greater role in the quantitative understand- 
ing of small-scale perturbations observed in the lower 
stratosphere and in the prediction of disturbances at 
altitudes above the present database for aircraft de- 
sign criteria. The challenges and recommendations for 
further progress in the definition of strong discrete per- 
turbations at stratospheric altitudes are clear. Mete- 
orological observations and in-situ atmospheric data 
from the troposphere and lower stratosphere are re- 
quired to provide the basic statistical characterization 
of perturbation activity at these levels. Definition of 
perturbations at higher altitudes requires numerical 
simulation of the upward propagation dynamics. Such 
definition would help in determining how these pertur- 
bations are affected by the atmospheric structure and 
under what meteorological conditions instabilities will 
result in strong perturbations and turbulence at higher 
altitudes. As numerical simulation studies progress, 
their realism can be evaluated in stages by compari- 
son with in-situ perturbation measurements, with wave 
amplitude observations by rawinsonde balloon rise rate 
in the middle stratosphere, and with remotely detected 
wave activity in the mesosphere. Such work should 
significantly improve the accomplishment of two basic 
goals. These goals are as follows: 

• To predict the nature of wave propagation and in- 
stabilities at higher altitudes 

• To improve estimates of the perturbation-length 
scales and maximum magnitudes for hypersonic 
aircraft design purposes 

These accomplishments will inherently support the def- 
inition of hypersonic vehicle design criteria and the for- 
mulation of simple, discrete models for characterizing 
perturbations in the vertical and horizontal wind com- 
ponents as well as the pressure, density, and tempera- 
ture states. 


Discrete Perturbation Model Concept 

This section describes a simple, initial, generic con- 
cept for a combined discrete atmospheric perturbation 
model. Such models may be used for design sensitiv- 
ity assessment. However, their use as design criteria 
throughout the stratosphere will not be warranted un- 
til the appropriate magnitudes, scale lengths, and risk 
levels are better estimated for the higher altitudes. 

The traditional derived equivalent gust velocity, 
is basically a discrete perturbation form which has 
served its purpose well as an airframe lifetime gust 
loads predictor for altitudes which were surveyed with 
similar aircraft. 5 Its one-minus-cosine load shape, in 
which the gust encounter builds up in the positive di- 
rection and decays back to zero over a distance of 25 
wing chord lengths, is a one-sided model for the verti- 
cal gust component alone. As airframe structures be- 
came more flexible, it became prudent to character- 
ize gust data in terms of the turbulence PSD for the 
evaluation of second-order responses. 4 For some pre- 
liminary design studies, simple, repeatable inputs of 
sharp, two-sided transients which represent the larger 
magnitude perturbations embedded in continuous tur- 
bulence are desirable. Hypersonic vehicle airworthi- 
ness may depend on gust-induced responses for several 
subsystems, such as propulsion, flight control, and ba- 
sic sensor capabilities, in addition to the gust loads 
on the structure. Therefore, the discrete perturbation 
model must include all physical components of the at- 
mospheric transients that may affect the vehicle. 

One simple, generic model that combines all of the 
perturbation variables is a vortex with solid body ro- 
tation in the core and velocity components outside 
the core decaying inversely with distance from the 
core. This perturbation form is generally characteris- 
tic of atmospheric observations for phenomena includ- 
ing vortex-ring segments accompanying downburst out- 
flows, dust devils, and higher altitude rotor formation 
resulting from Kelvin- Helmholtz instability and other 
familiar vortex shedding phenomena . 21,47 Generation 
of the instability cells results from wave amplification 
which brings up colder, more dense air and brings down 
warmer, less dense air. Figure 13 shows potential tem- 
perature isotherm distortion leading to formation of a 
strong turbulent perturbation. 

When vortex formation occurs with the axis of the 
core in the horizontal plane, the gust velocity com- 
ponent perpendicular to the flightpath achieves maxi- 
mum perturbation velocity along trajectories through 
the center of the core. The perturbation component 
parallel to the flightpath attains maximum amplitude 
for trajectories tangential to the outer edge of the core. 
For illustration of the resulting perturbations, the max- 
imum gust velocity of 10 m/sec with the vortex core 
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radius dimensioned to 200 m is selected. This com- 
bination of maximum gust velocity and core radius 
produces an arbitrarily strong maximum shear of 0.05 
sec" 1 . 

Figure 14(a) shows the vertical gust velocity for a 
horizontal flight trajectory through the center of the 
core. In comparison, note that the peak- to- peak per- 
turbation amplitude through the core is approximately 
the same as prescribed for maximum cruise speed at 
15 km. 1 The horizontal component for this trajectory 
directly through core center of the idealized pertur- 
bation does not experience any variation. Both hor- 
izontal and vertical components experience compara- 
ble perturbation magnitudes for horizontal trajectories 
through the upper or lower edges of the core. 

Figure 14(b) illustrates the horizontal and vertical 
components of the discrete perturbation model for a 
trajectory tangential to the (that is, through the) up- 
per edge of the core. Again, note that the amplitudes 
of this discrete model are approximately the same as 
the Ref. 1 criteria for design loads criteria at maxi- 
mum dive speed at an altitude of 15 krn. For the tem- 
perature change, a value of 10 K with buildup and 
drop-off distances of 400 m is used. This value ap- 
pears to be reasonably conservative on the basis of 
previously reported aircraft measurements. 48 Because 
temperature transient peaks often do not coincide with 
peaks in the gust velocity components, the coldest tem- 
perature is ofTset 200 m to the updraft side of the core 
as illustrated with the attending ambient density vari- 
ation in Fig. 14(c). This configuration of cold air rel- 
ative to the motion components is highly arbitrary be- 
cause neither the most prevalent nor the most critical 
patterns are known. Pressure is approximated by use 
of Bernoulli’s equation and the equilibrium assumption 
that the pressure gradient force balances the centrifu- 
gal force. For this model, the pressure perturbation 
is less than 0.1 percent at the outer edge of the core. 
This perturbation is not illustrated because it results 
in a nearly negligible effect on density. 

Specification of such simple, discrete perturbation 
models as well as their geometric dimensions and per- 
turbation magnitudes is reasonably straightforward for 
the altitudes of available in-situ data below approxi- 
mately 20 km. At higher altitudes above available in- 
situ data, selection of the gust and thermodynamic per- 
turbation magnitudes as well as the appropriate geo- 
metric length scales becomes highly questionable. Dis- 
crete perturbation model specification for hypersonic 
aircraft design at these higher altitudes requires the 
use of numerical simulation of the small-scale dynamic 
atmospheric processes which lead to strong perturba- 
tions in gust velocities, temperature, and density. 


Concluding Remarks 

Vehicle design criteria for lower altitudes have served 
their purpose well for structural design of subsonic air- 
craft to altitudes of approximately 15 km. However, 
these criteria do not incorporate combined perturba- 
tion inputs for either the motion components or the 
temperature, pressure, and density states. These states 
are expected to be more significant to airbreathing 
hypersonic vehicles. In addition, the amount of in- 
situ data for establishing gust design criteria decreases 
markedly between altitudes of 15 and 25 km. Above 
these altitudes, essentially no in-situ data have been 
gathered. Available data for turbulence and temper- 
ature transients encountered by aircraft in the lower 
stratosphere often show an association with lower alti- 
tude mountain-wave activity. Therefore, improved un- 
derstanding of upward wave propagation processes is a 
key element in the formulation of atmospheric pertur- 
bation design criteria for higher altitude aircraft. 

Recent advances in numerical simulation of wave 
propagation processes are making such improvements 
in understanding the conditions and atmospheric struc- 
tures associated with the development of wave in- 
stabilities which cause turbulence and strong gusts. 
A current experiment comparing in-situ observa- 
tions of mountain-wave-induced perturbations in the 
lower stratosphere with two-dimensional, numerical, 
atmospheric dynamics simulations made by separate 
codes was described. Initial results from these com- 
parisons indicate numerical simulations will provide 
useTul descriptions of higher altitude perturbations. 
An example of a simple, generic, discrete perturba- 
tion model combining motion components and ther- 
modynamic (pressure, density, and temperature) dis- 
turbances was given. Definition of appropriate distur- 
bance magnitudes and length (dimensional) scales for 
such discrete atmospheric perturbation models at alti- 
tudes of 25 to 50 km is needed for hypersonic vehicle de- 
sign. These perturbation magnitudes and length scales 
can be defined by studies which combine these numer- 
ical atmospheric simulation tools with higher altitude 
observations in the middle and upper stratosphere. 
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Fig. 1 Representative atmospheric structure, January average for Edwards, California (Ref. 15). 
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Fig. 2 High-altitude turbulence survey results for amount of turbulence expressed as the percent of flight distance 
(Ref. 13). 
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Fig. 3 Cumulative frequency distributions of peak accelerations caused by turbulence. 
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Fig. 4 Seasonal variation of turbulence encountered in the lower stratosphere between 12 and 17 km. 
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(a) Rate variations of 50 to 150 m/min for turbulence (b) Rate variations of 50 to 150 m/min for turbulence 
near Ely, Nevada, near Winnemucca, Nevada. 
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(c) Rate variations of 30 to 70 rn/min for nonturbu- 
lence near Las Vegas, Nevada. 
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(d) Rate variations of 30 to 70 m/min for nontur- 
bulence near Ely, Nevada. 


Fig. 5 Rawinsonde balloon ascent rate profiles for turbulence and nonturbulence areas encountered in high-altitude 
supersonic flight (Ref. 18). 
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Fig. 6 Mountain- wave turbulence forecast verification (modified from Ref. 19). 
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Fig. 7 Cross-section of potential temperature and wi 
turbulence (Ref. 28). 
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Fig. 8 Input temperature profile and rawinsonde observation data for simulation of mountain- wave turbulence in 
the stratosphere. 
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Fig. 9 Input wind profile and rawinsonde observation data for mountain-wave turbulence in the stratosphere. 
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(a) Potential temperature contours (K, deg) for strong turbulence. 



(b) Horizontal velocity contours (m/sec) for strong turbulence. 

Fig. 10 Vertical and horizontal cross-sections for numerical simulations of potential temperature and wind contours 
with and without observed turbulence in the stratosphere. 
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(c) Potential temperature contours (K, deg) without turbulence. 



(d) Horizontal velocity contours (m/sec) without turbulence. 
Fig. 10 Concluded. 
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(c) Well-developed turbulence regions, 14 hr. (d) Waves and turbulence decay, 16 hr. 

Fig. 11 Numerical simulation of wave amplitude growth, breakdown and overturning, instability, and decay in the upper stratosphere (Ref. 43). 
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(a) At 4 hr, monotonic wave upward propagation approaches the critical level at 25 km. 
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(b) At 9 hr, critical level interaction induces zones of reverse flow between 22.5 and 25 km. 

Fig. 12 Numerical simulation stream function field showing instability development beneath the critical level at 
25 km resulting from upward propagation of a monotonic wave at an altitude of 15 km (Ref. 44). 
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(c) At 10 hr, smaller scale perturbations develop. 
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(d) At 12 hr, perturbations have developed at smaller scales and propagate to both higher and lower altitudes. 


Fig. 12 Concluded. 
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Fig. 13 Schematic of distorted atmospheric potential temperature stratification leading to instability and vortex 
wrap up. 
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(a) Vertical component for horizontal trajectory through the core center. 



(b) Components for a trajectory through the core top. 



(c) Thermodynamic perturbations for a trajectory through the core top. 

Fig. 14 Discrete, small-scale, two-dimensional, atmospheric perturbation model. 
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